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I.  INTRODUCTION 


The  objective  of  the  research  described  In  this  report  was 
to  investigate  the  potential  use  of  scattered  ultraviolet  radi- 
ation In  the  solar  blind  region  of  the  spectrum,  as  a means  of 
voice  communication.  A communication  system  using  scattered  UV 
light  would  utilize  a frequency  modulated,  high  repetition  rate, 

UV  flash  lamp  as  a transmitter  and  a combination  of  solar  blind 
filters  and  photomultiplier  tube  as  a receiver.  The  voice  com- 
munication could  thus  be  accomplished  day  or  night,  bad  weather 
or  good,  over  hills,  around  corners  or  line  of  sight.  The  system 
would  be  inherently  short  range  (due  to  ozone  absorption),  covert, 
secure.  Jam  resistant.  Invisible  to  the  eye  and  would  not  require 
a conventional  antenna. 

Virtually  no  sunlight  reaches  the  earth's  surface  in  the  solar 
blind  region  of  the  spectrum  (2200  to  2800^)  due  to  ozone  absorp- 
tion in  the  upper  atmosphere.  Consequently,  an  optical  communica- 
tion system  can  operate  without  interference  from  the  sun  in  this 
spectral  region.  Typical  distributions  of  ozone  concentration  as 
a function  of  altitude  are  given  in  Figure  1.  It  car,  be  seen  that 
the  ozone  concentration  peaks  at  approximately  25  km  and  generally 
has  a rather  low  value  at  sea  level.  The  relatively  low  concentra- 
tion of  ozone  at  the  earth's  surface  permits  transmission  of  ultra- 
violet radiation  over  reasonable  distances  and,  at  the  same  time, 
has  Important  implications  with  respect  to  security  and  antijamming 
characteristics . 

Transmitted  ultraviolet  radiation  undergoes  strong  Mle  type 
scattering  at  low  altitudes  due  to  the  high  scattering  cross-sections 
of  atmospheric  aerosols  in  this  wavelength  region.  During  periods 
of  bad  weather  the  concentration  of  atmospheric  aerosols  is  likely 
to  Increase,  thus  enhancing  the  Intensity  of  scattered  radiation. 
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Figure  1.  Variation  of  ozone  density  distribution  with  altitude, 
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It  3s  through  the  scattering  mechanism  that  transmitted  UV  light 
Is  able  to  travel  over  and  around  natural  and  manmade  obstacles. 

The  nature  of  the  scattering ' process  In  the  ultraviolet  Is  therefore 
very  important,  primarily  with  respect  to  preferential  scattering 
In  the  forward  and  backward  direction.  Unfortunately,  virtually 
no  fundamental  scattering  data  are  available  for  the  solar  blind 
region  of  the  spectrum.  An  experimental  and  theoretical  investi- 
gation of  the  manner  in  which  the  atmosphere  processes  a single 
flash  of  ultraviolet  radiation  from  an  isotropic  source  has  there- 
fore been  undertaken  as  one  of  the  major  objectives  of  this  study. 

The  individual  goals  of  this  investigation  described  in  this 
report  can  be  summarized  as  follows: 

1.  To  determine  the  optimum  components  currently 
available  for  an  ultraviolet  communication 
system  Including  solar  blind  filters,  photo- 
multiplier tubes,  high  frequency  UV  flash- 
lamps,  triggering  circuits  and  other  associated 
electronics . 

2.  To  design  ana  fabricate  a proof-of-prlnciple 
transmitter  and  receiver,  over  which  voice 
communication  can  actually  be  accomplished  via 
scattered  ultraviolet  radiation. 

3-  To  test  the  efficiency  of  UV  flashlamps  In  con- 
verting input  power  into  UV  energy,  and  to  test 
the  effects  of  various  filter  combinations  and 
field  of  view  ffects  on  the  received  signal. 

4.  To  determine  experimentally  the  manner  in  which 
the  atmosphere  processes  (scatters  and  absorbs) 

UV  radiation. 
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5.  To  mathematically  model  the  observed  atmospheric 
scattering  phenomena. 

6.  To  combine  the  knowledge  obtained  from  all  of 
the  above  to  make  an  estimate  of  the  ultimate 
performance  and  power  requirements  of  a practical 
UV  communications  sytem. 

7.  To  outline  the  advantages  and  disadvantages  of  a 
UV  communication  system  and  to  envision  potential 
applications . 

Finally,  it  is  the  purpose  of  this  effort  to  Identify  the  areas 
where  furth.er  research  needs  to  be  carried  out  before  a practical 
system  can  be  developed. 
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II.  APPARATUS 


Ultraviolet  Sources 

Two  types  of  ultraviolet  sources  were  employed  In  the  present 
study.  Small  bulbs  with  high  pulse  repetition  rates  were  used  in 
the  Investigation  of  voice  communication.  Large,  high  power  flash- 
lamps  were  used  for  obtaining  measurements  of  atmospheric  scatter- 
ing. Both  types  of  lamps  generate  ultraviolet  radiation  by  pass- 
ing electric  arcs  through  hydrogen-xenon  mixtures.  Typical  spectral 
distributions  for  these  lamps  are  shown  In  Figure  3- 

The  high  power  flashlamps  were  obtained  from  the  EG4G  Company. 
These  bulbs  (FX-36C-3)  have  a maximum  output  power  of  1000  Joule£ 
per  flash.  A photograph  of  one  of  these  lamps  mounted  on  the 
triggering  unit  Is  shown  In  Figure  ^a.  Figure  ^<a  does  not  show 
the  large  set  of  capacitors  which  are  employed  to  provide  the  pulse 
energy  to  the  flashlamp.  The  set  contains  three  capacitors  rated 
at  80  uF  - 2000  V-DC  and  one  capacitor  rated  at  200  uF  - 2500  V-DC 
for  a total  rating  of  ^^0  pF. 

2 

Since  the  power  of  the  flashlamps  Is  0.5  CV  where  C is  the 
capacitance  and  V Is  the  voltage,  the  power  may  be  doubled  by  in- 
creasing the  voltage  by  /5”  . Thus,  voltages  of  1^14,  1000,  707 » 

and  500  hav  been  used  for  the  atmospheric  scattering  studies.  The 
basic  electrical  circuit  for  the  flashlamp  is  shown  In  Figure  5* 

With  this  circuit  the  pulsewldth  of  the  lamp,  or  the  duration  of 
the  flash.  Is  proportional  to  / LC  where  L Is  the  circuit  flnduct- 
ance.  Typical  pulse  profiles  are  shown  In  Figure  6.  For  the  440  pF 
set  of  capacitors,  the  flashlamp  produced  a pulse  of  approximately 
1 msec  duration. 

The  high  repetition  rate  bulbs  were  obtained  from  United  States 
Scientific  Instruments  Inc.  These  bulbs  were  of  two  basic  designs, 
one  rated  at  10  watts,  the  other  at  30  watts.  Both  bulbs  were  suc- 
cessfully operated  at  8000  cycles  per  second.  The  high  frequency 
trigger  circuit  for  these  bulbs  Is  shown  in  Figure  7»  The  10  watt 
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Figure  4. 


Photographs  of  the  UV  sources  employed  for  the 
scattering  measurements  (4a)  and  for  the  voice 
communication  experiments  (4b). 
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Figure  6.  Typical  pulst-;  profiles  for  high  eaiergy  flashlamp. 
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Figure  7.  High  frequency  trigger  circuit  for  conununication  source. 


bulb  Is  shown  In  Figure  4b  mounted  within  a box  containing  the 
trigger  circuit.  The  front  surface  of  the  lamp  Is  made  of  Corning 
type  7-54  glass  which  does  not  transmit  In  the  visible  spectral 
region.  The  30  watt  bulb  which  Is  not  shown  has  been  mounted  on 
top  of  the  box  shown  In  Figure  4b. 

An  experimental  analysis  of  the  two  bulbs  has  been  conducted 
to’ determine  the  relative  UV  output  as  a function  of  Input  power. 

It  should  be  indicated  at  this  point  that  the  bulbs  were  designed 
for  photocopying.  Thus,  the  ITV  radiation  has  been  Intentionally 
suppressed.  The  results  of  our  analysis  are  shown  In  Figures  8 and 
9.  As  can  be  seen  from  the  figux’es,  enhancement  of  the  U'V  radi- 
ation requires  high  voltage  and  low  capacitance  for  a given  power 
setting.  The  large  Increase  in  UV  output  with  power  is  shown  in 
Figure  9.  As  the  input  power  Is  Increased  from  4 to  6 watts,,  the 
UV  output  Increases  by  almost  300?.  Unfortunately,  the  lamps  and 
associated  trigger  circuitry  were  originally  designed  to  achieve 
the  rated  power  with  high  capacitance  and  low  voltage  which  enhances 
the  visible  radiation  used  in  photocopying.  At  voltages  above  190V 
shorting  occurs  within  the  lamp  circuitry  thereby  establishing  an 
upper  voltage  limit.  The  power  obtainable  from  A. R. A. P.'s  50  mllll- 
amp  DC  power  supply  Is  consequently  limited  to  less  than  10  watts  of 
average  power.  A conversion  efficiency  of  Input  power  to  UV  output 
of  approximately  .1?  has  been  measured  at-  the  highest  power  setting 
shown  In  Figure  9- 

Filters 

The  type  of  filters  necessary  to  provide  the  required  long  wave- 
length solar  radiation  rejection  are  quite  unique.  Historically, 
the  filters  have  been  composed  of  numerous  substances  such  as  nickel 
sulphate  hexahydrate,  catlon-x.  Corning  glass  9863,  and  aluminum 
Interference  filters.  The  various  transmlttances  of  these  materials 
and  several  combinations  are  shown  In  Figure  10.  The  basic  problems 
with  these  filters  are  the  low  peak  transmittance,  3-5^,  and  the 
size,  minimum  1.5-lnches  total  depth.  These  filters  currently  are 
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Figure  8.  Relative  energy  conversion  of  USSI  bulbs. 
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Power  response  of  30  watt  USSI  bulb. 


Figure  9 . 
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Figure  10.  Transmittance  curves  of  materials  emp].oyed 
in  solar  blind  systems. 
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provided  by  Barr  Associates,  Cambridge,  Massachusetts.  Filters 
have  been  developed  by  the  author  (Ref.  2),  in  collaboration  with 
Optical  Coating  Laboratories,  Inc.  which  were  0.25-inches  thick 
but  with  a peak  transmittance  of  1.555  at  2630X.  Under  the  present 
contract,  five  filter  combinations  have  been  obtained  from  Barr 
Associates  to  provide  considerable  flexibility  in  selecting  several 
wavelength  intervals  between  220oX  and  2800X  with  considerable  vari- 
ations in  long  wavelength  rejection.  The  transmission  chr'^acter- 
Istlcs  of  these  filters  are  indicated  in  Figures  11  and  12.  Filters 
1 and  2 are  the  basic  filters  necessary  to  obtain  a high  degree  of 
solar  blindness.  The  other  three  filters,  used  with  the  basic 
filters,  provide  wavelength  discrimination  in  the  spectral  region 
between  2^40oX  and  2700X.  When  these  supplemental  filters  are  com- 
bined with  the  PMT  and  the  basic  filters,  we  obtain  the  overall 
system  spectral  response  shown  in  Figure  13.  These  elements  are 
shown  in  Figure  l4 . 

An  additional  0.5-lnch  piece  of  NlSOj^6(H20)  has  been  used  in 
conjunction  with  the  basic  filters  nos.  1 and  2 and  the  2680  filter. 
The  purpose  of  this  additional  filter  is  to  eliminate  the  effects  of 

O 

a light  leak  at  4000A  associated  with  an  "off  design"  photomulti- 
plier tube  described  in  the  next  section.  Although  the  NlBOj^SCH^O) 
filter  has  a transmittance  of  approximately  95!f,  the  net  effect 
Including  a reduction  in  f leld-of-view  is  a 2155  decrease  in  signal. 
The  field-of-vlew  function  (FOV)  for  this  combination  of  filters 
has  been  determined  experimentally  by  observing  a point  source  of 
light  at  various  angles  from  the  detector  optical  axis.  The  results 
are  shown  in  Figure  15  and  have  been  Included  in  the  simplified 
scattering  model  described  in  a later  section.  The  total  fleld-of- 
view  is  seen  to  be  approximately  28°  with  a half-angle  (defined 
as  the  half-angle  at  50!(  maximum  intensity)  of  7®. 

For  the  purpose  of  measuring  the  atmospheric  attenuation  co- 
efficient, it  is  necessary  to  point  the  detector  directly  at  the 
source  with  a /ery  narrow  fleld-of-vlew.  The  purpose  of  the  narrow 
fleld-of-vlew  is  to  separate  transmitted  light  from  scattered  light. 
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Figure  11.  Transmission  characteristics  of  basic 
filters  1 and  2. 
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Figure  1^4.  Photograph  of  components  of  filter  assembly. 


Figure  15.  Angular  dependence  of  field  of  view  function. 
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A quartz  lens  with  a four-inch  focal  length  was  therefore  attached 
to  the  front  of  the  filter  assembly.  Four  Inches  behind  the  len;  , 
a black  aluminum  disk  with  a l/l6-lnch  diameter  hole  was  Inserted. 

This  combination  of  lens  and  stop  yield  a flat  field  of  view 
function  with  a total  angle  of  only  1.8°.  For  atmospheric  scatter- 
ing measurements,  the  stop  and  lens  were  removed. 

Photomultiplier 

The  basic  photomultiplier  employed  In  this  Investigation  was 
the  5^2P-09-18  provided  by  EMR.  This  tube  has  a rubidium  tellurlde 
photosensitive  cathode  with  the  quantum  efficiency  as  given  In 
Figure  16.  The  particular  advantage  v)f  this  tube  Is  the  very  low 
dark  count  (the  random  emission  of  electrons  by  the  photosensitive 
cathode)  which  Is  approximately  10  counts/sec.  A very  low  dark 
count  Is  Important  In  the  measurement  of  atmospheric  scatteilng  since 
we  can  measure  much  lower  levels  of  radiation.  In  addition,  since 
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this  tube  provides  a very  large  dynamic  range,  10  , very  high  levels 
of  radiation  can  be  measured. 

Detection  System 

The  basic  electronic  circuitry  has  been  designed  around  the 
detection  of  the  weakest  signals  from  the  ultraviolet  source.  The 
pulse  of  electrons  produced  at  the  anode  of  the  photomultiplier, 
which  is  Initiated  by  the  ejection  of  an  electron  from  the  photo- 
sensitive cathode,  typically  has  a pulse  width  of  ^0  nanoseconds  (ns). 
This  type  of  pulse  Is  ideally  suited  for  digital  logic  data  processing. 

The  basic  detection  system  Is  shown  in  Figure  17  which  also 
indicates  the  various  components.  Figure  17a  Is  a photograph  of 
the  detector  and  the  HP  counter  employed  for  photon  counting.  The 
filter  assembly  also  Is  shown.  Figure  17b  Indicates  the  various  com- 
ponents comprising  the  detector.  The  EMR  photomultiplier  (PMT)  is 
a replacement  for  the  original  PMT  which  rapidly  lost  sensitivity 
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Wavelength,  nanometers 


Figure  16.  Measured  quantum  efficiency  for  EMR  5^2P-09-l8 
photomultiplier , 


2) 


Figure  17*  Photographs  of  the  basic  detecting  system 
indicating  pertinent  components. 


over  a period  of  one  day.  The  replacement  PMT  does  not  have  the 
long  wavelength  rejection  normally  associated  with  EMR  rubidium 
teiluride  PMT's.  However,  EMR  did  not  plan  to  make  any  new  tubes 
in  the  near  future.  Therefore,  they  searched  their  archives  con- 
taining tubes  which  did  not  meet  their  very  tight  requirements. 

The  present  tube  evidently  did  not  meet  their  required  long  wave- 
length rejection.  Although  some  radiation  around  4000a  can  be 
detected,  the  amount  is  small  and  has  not  hindered  our  measure- 
ment effort.  The  EMR  power  supply,  which  converts  28V-DC  to  2000- 
36OOV-DC,  also  failed  but  was  replaced  by  EMR. 

The  pulse  amplitude  discriminator,  PAD,  provides  a very  large 
range  in  sensitivity.  This  PAD  amplifies  the  40ns  (nanoseconds) 
pulses  from  the  PMT  and  produces  a positive  pulse  of  60ns  and 
3.5  V-DC.  The  maximum  random  count  rate  is  4x10^  cts/second. 

Since  the  basic  mode  of  measurement  is  that  of  pulse  or  photon 
counting,  all  measurements  indicate  the  number  of  photons/sec 
detected  by  the  system.  This  type  of  measurement  permits  a wide 
dynamic  range.  All  measurements  are  in  terms  of  cts/sec  or  cts/tlme 
Interval.  The  dark  count  of  the  system  is  the  number  of  counts 
received  per  second  with  the  entrance  aperture  completely  closed. 

The  count  can  be  varied  somewhat  by  the  sensitivity  level  set  by 
the  PAD.  Under  normal  conditions  the  dark  count  ranges  between 
6 and  15  cts/sec. 

The  apparatus  which  was  used  to  obtain  field  measurements  of 
scattered  radiation  is  shown  in  Figure  I8.  The  table  mounted  on 
a tripod  contains  the  detection  system,  a telescope  for  sighting 
on  the  source,  a compass  used  to  obtain  bearing  and  azimuthal 
angle  from  source,  and  a device  to  measure  elevation  angle. 

Transmitter  and  Receiver 

Voice  communication  has  beeri  achieved  employing  a breadboard 
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transmitter  and  receiver.  The  system  employs  pulse  frequency 
modulation  with  a center  frequency  of  4810  pps.  The  voice  band- 
width is  restricted  to  the  Interval  200  Hz  - 2000  Hz.  Figures 
19  and  20  indicate  block  diagrams  of  the  associated  electronic 
circuits . 

This  particular  system  Includes  a feature  which  helps  to 
eliminate  random  pulses.  This  is  achieved  by  a minimum  Interval 
detector.  After  receiving  a pulse  the  system  will  not  admit 
another  pulse  of  l67psec.  Then  the  system  is  "on"  for  83ysec  to 
admit  the  next  pulse.  In  this  manner,  extraneous  pulses  tend  to 
be  eliminated. 

At  the  present  time,  the  transmitter  circuit  does  not  have 
all  of  the  elements  indicated  In  Figure  19*  The  automatic  gain 
control,  AGC,  and  the  logarithmic  speech  compressor  have  not 
been  Installed.  The  receiver  circuit  has  all  of  the  indicated 
elements . 

Our  basic  detection  system,  Including  the  EMR-PMT,  was  used 
as  the  front  end  of  the  receiver.  The  60  ns  output  pulses  were 
too  sharp  for  the  input  buffer  circuit.  Therefore,  a fast  one- 
shot  multivibrator  was  used  to  stretch  the  pulse  to  6ysec. 
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III.  MEASUREMENTS 


The  purpose  of  the  field  measurements  was  to  evalua.  tl.u 
potentia.1  of  nonline  of  sight  communication  via  scattered  ulti’a- 
violet  radiation  in  the  lower  atmosphere.  A primary  concern  in 
the  evaluation  of  the  communication  system  Is  the  possible  inter- 
ference which  might  arise  due  to  natural  and  manmade  UV  sources. 

Thus,  the  first  phase  was  the  determination  of  the  ambient  back- 
ground. The  next  step  was  the  determination  of  the  manner  In 
which  the  atmosphere  processes  a flash  of  UV  photons.  The  magni- 
tude of  UV  scattering  and  absorption  by  the  atmosphere  was  deter- 
mined by  surveying  the  sky  from  various  locations  away  from  the 
source.  In  addition,  to  correlate  the  data  with  our  mathematical 
modeling,  it  was  necessary  to  determine  experimentally  what  portion 
of  the  received  signal  came  from  the  lower  part  of  the  field  of 
view  and  what  portion  came  from  the  upper  part.  Finally,  field 
measurements  provided  first-hand  experience  In  determining  problems 
which  might  arise  with  an  actual  UV  communications  system. 

Background  Measurements 

The  most  significant  background  measurements  were  those  obtain- 
ed during  relatively  clear,  sunny  days  In  which  solar,  radiation  was 
expected  to  be  the  dominant  source  of  scattered  UV  radiation.  Numer- 
ous measurements  were  made  employing  various  filter  combinations. 
These  measurements  Included  those  in  which  the  sun  was  directly  in 
the  detector  field  of  view.  The  results  are  indicated  in  Table  1. 
This  table  contains  two  sets  of  measurements  approximately  one  hour 
apart.  The  bearing  and  elevation  of  the  sun  corresponded  to  the 
beginning  of  the  first  set  of  measurements.  Because  of  the  long 
wavelength  "leak"  around  400oX,we  used  a long  wavelength  filter, 

O 

LWP,  which  was  opaque  to  radiation  below  approximately  3700A.  To 
determine  the  UV  radiation  detected  by  the  system,  we  would  make 
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TABLE  1 

TYPICAL  RESULTS  PROM  SKY  SURVEY 

SUN  AT  BEARING  255°  ELEVATION  48° 

SYSTEM  DARK  COUNT  5.2  CTS/SEC  OVER  50  SECONDS 


FILTER  [ 

COMBINATION  BEARING 

ELEVATION 

AVERAGE 

CTS/SEC 

TIME  PERIOD 
SEC 

uv 

CTS/SEC 

r 

1,2 

45 

20 

44.7 

50 

31.3 

l,2,LWx» 

45 

20 

13.4 

30 

1,2 

90 

20 

41.0 

40 

25.3 

1,2,LWP 

90 

20 

15.7 

30 

1,2 

225 

20 

55.1 

30 

32.9 

1,2,LWP 

225 

20 

22.2 

30 

1,2 

315 

20 

63.9 

30 

44.0 

1,2,LWF 

315 

20 

19.9 

20 

1,2 

270 

20 

76.1 

40 

47.3 

1,2,LWF 

270 

20 

28.8 

20 

1,2,2680 

270 

20 

i4.3 

30 

0 

1,2, 2680, LWF 

270 

20 

i4.5 

30 

1 

1,2 

DIRECTLY  AT 

SUN 

52.3 

70 

21.1 

1,2, LWF 

11 

1! 

ft 

31.2 

50 

1,2,2680 

11 

If 

11 

i4.5 

60 

0 

1,2,2580 

11 

11 

11 

16.8 

90 

0 

1,2,2480 

11 

11 

11 

17.2 

40 

0 

I 1,2, 2480, LWF 

It 

11 

11 

15.7 

60 

0 

»‘LWP  is  a filter  transparent  at  wave 


lengths  longer  than  3700X. 


31 


one  set  of  measurements  without  the  LWF  and  on<  set  with  the  LWP. 
Thus,  the  detected  UV  radiation  would  be  the  difference  between 
the  two  measurements.  This  value  Is  In  the  last  column.  Filters 

O 

1 and  2 do  permit  some  solar  radiation  lii  the  28C0-2900A 

spectral  region  to  reach  thj  PMT.  This  problem  was  anticipated 

and,  In  fact,  the  filters  were  designed  in  that  manner.  However, 

o o o 

when  any  of  the  2680A,  258OA,  and  2480A  filters  were  used  with  1 
and  2, no  solar  UV  radiation  was  detected.  The  counts  received 
when  these  filters  were  employed  corresponded  to  the  ''leak"  around 

O 

4000A.  The  last  series  of  measurements  with  1,  2,  and  2480  with 
and  without  the  LWB'  clearly  indicates  that  no  solar  UV  radiation  is 
being  detected.  The  count  level  of  approximately  14-17  is  composed 
of  5.2  counts  of  average  dark  count  and  the  remainder  from, the 

O 

leak  at  4000A.  In  all  measurements  after  this  Initial  set,  one  of 

00  O 

the  three  filters  2680A,  258OA,  and  2480A  was  slways  employed 
with  1 and  2. 

O 

The  basic  filter  combination  of  1,  2,  and  2680A  was  extensively 
employed  In  measurements  of  tae  natural  and  manmade  backgrounds  In 
and  around  Princeton  ana  Tre  .ton,  N.J.  Considerable  time  was 
devoted  to  this  task.  Surveys  were  made  during  daylight  hours  and 
early  evening  hours.  The  early  evening  hours  Included  times  during 
which  street  lights  and  neon  signs  were  on.  The  count  rate  never 
at  any  time  Indicated  a source  of  UV  radiation  In  the  spectral 

O 

region  below  280OA.  Counts  were  registered  due  to  noisy  transformers 
In  the  necn  lights  but  no  UV  radiation  was  detected.  This  result 
was  not  surprising  since  all  manmade  sources  of  Illumination  have 

O 

glass  covers  which  do  not  radiate  below  2800A. 
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Measurements  have  been  made  of  possible  UV  radiatioi  emitted 
by  lightning  during  thunderstorms.  A number  of  lightning  strokes 
were  observed  with  ranges  from  0,6  km  to  km.  Low  counts  were 
observed,  In  the  400-600  count  region,  but  these  were  most  probably 
due  to  the  electrical  disturbance  produced  by  the  lightning.  Simi- 
lar counts  were  recorded  with  the  aperture  open  and  closed.  Thus, 
the  UV  radiation  which  may  be  produced  by  lightning  strokes  Is 
negligible . 

As  expected,  the  ultraviolet  radiation  emitted  from  flames 
was  detectable  with  the  UV  system.  The  flames  from  matches,  lighters 
and  propane  torches  were  used  on  several  occasions  to  test  the 
vulnerability  of  the  UV  communication  system  to  background  noise. 
Unfortunately,  no  opportunity  arose  during  the  course  of  the  measure- 
ments to  observe  a large  fire  such  as  a burning  building.  Nor  did 
an  opportunity  arise  to  measure  the  effects  of  explosives  and 
munitions . 

High  Energy  Flas hlamp  Measurements 

Two  basic  types  of  measurements  were  made  using  the  high  energy 
flashlamp.  By  pointing  the  detector  directly  at  the  source  from 
various  distances,  the  llne-of- sight  attenuation  coefficient  was 
measured.  The  attenuation  coefficient  Is  a direct  measure  of  the 
number  density  and  cross-sections  of  atmospheric  absorbers  and 
scatterers  on  a given  dry.  Alternately,  the  detector  was  pointed 
away  from  the  UV  source,  so  as  to  systematically  determine  the 
amount  of  scattered  radiation  coming  from  various  portions  of  the 
sky.  These  maps,  which  indicate  the  way  In  which  UV  radiation  Is 
processed  by  the  atmosphere,  are  referred  to  as  scattering  "foot- 
prints." 

Attenuation  Measurements 

In  principle,  the  measurement  of  attenuation  coefficients  Is 
straightforward.  The  attenuation  coefficient  defined 
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as  the  fraction  of  light  removed  or  deflected  from  Its  original 
direction  due  to  absorptloii  and  scattering,  per  unit  path  length 
at  wavelength  X , In  the  case  of  radiation  passing  through  a 
uniform  atmosphere  from  an  Isotropic  source,  the  relation  between 
the  transmitted  and  Initial  Intensity  Is 


I(R) 


(1) 


where  the  Initial  source  Intensity  Iq  Is  In  units  of  photons  per 
flash,  I(R)  Is  In  units  of  photons/km^/flash,  and  Is  In  units 

of  1/km.  If  0^  Is  equal  to  zero,  then  the  above  relation  reduces 
to  the  simple  Inverse  square  law.  If  the  area  of  the  detector 
entrance  optics  Is  and  the  conversion  efficiency  of  captured 

photons  to  photon  counts  Is  POV  (the  fleld-of-vlew  function),  then 
the  photon  count  PC  due  to  a single  flash  of  photons  at  wavelength 
(X)  is 


PC(X,R) 


(FOV)ApI  ^ 

HttR^ 


(2) 


Because  an  accurate  value  for  the  combined  term  (FOV)Aj^I^  Is  not 
always  known,  the  value  of  usually  obtained  by  measuring 

PC(X)  at  two  different  distances  R^  and  R2  • The  attenuation  co- 
efficient can  then  be  solved  from  the  relation 

, PC(R  X) 

" spiq;  ^ pftdij, n 

In  the  present  study,  measurements  have  not  been  made  at  discrete 
wavelengths  but  rather  over  narrow  bands  such  as  those  shown  In 
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Figure  13.  Therefore,  It  has  been  assumed  that  the  attenuation 
coefficients  measured  using  the  above  equation  are  average  values 
over  the  width  of  the  band. 

The  attenuation  coefficient  actually  represents  both  the 
decrease  in  total  UV  radiation  due  to  ozone  absorption  and  the 
deflection  of  UV  radiation  due  to  scattering.  The  basic  dif- 
ference between  absorption  and  scattering  is  that  with  scattering 
the  photons  are  not  removed  from  the  atmosphere  but  merely  change 
direction.  The  attenuation  coefficient  can  be  written  as  the 

sum  of  the  absorption  coefficient  and  the  scattering  coefficient. 

°A  °ABS  ■*'  ^SCATT 

It  is  important  to  the  successful  measurement  of  the  attenua- 
tion coefficient  that  scattered  light  be  prevented  from  entering 
the  detector  when  the  light  source  is  viewed  head  on.  For  this 
reason,  a very  narrow  detector  field  of  view  is  required,  because 
scattered  light  will  reach  the  detecto  from  all  angles  while 
transmitted  light  comes  virtually  from  a "single"  angle  (see  Figure 
21).  While  it  is  not  possible  to  completely  eliminate  scattered 
light,  it  is  possible  to  reduce  it  to  a level  where  it  will  be 
small  compared  to  the  transmitted  light. 

Actual  values  of  the  attenuation  coefficient  as  measured  in 
the  field  varied  from  1.6/km  to  3* 2/km  which  is  consistent  with 
the  values  reported  by  Baum  and  Dunklernan  (Ref.  3)  as  shown  in 
Figure  22.  Early  in  the  study,  however,  the  values  of  the  attenua- 
tion coefficients  obtained  between  the  various  sites  on  a given 
day  were  not  consistent.  An  investigation  into  the  origins  of  the 
Inconsistency  showed  that  errors  in  the  attenuation  coefficient 
could  result  from  the  following: 
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Figure  21.  Reduction  of  scattered  light  for  measurement  of  attenuation  coefficient 
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Figure  22.  Baum  and  Dunkelman  attenuation  data. 
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incorrect  distance  measurements 


changing  visibility  during  the  course  of  the 
measurements 

imprecise  pointing  of  the  detector  toward 
the  source 

interference  from  scattered  radiation 
changing  flashlamp  intensities 
scintillation  error 

change  of  attenuation  (Coefficient  with  distance 

The  results  of  the  analysis  showed  that  while  interference  from 
scattered  radiation  can  result  in  the  measured  attenuation  co- 
efficient being  lower  than  the  actual  value,  the  low  values  should 
be  consistent  with  each  other.  Although  the  measured  coefficients 
can  be  very  sensitive  to  errors  in  the  distance  measurements  and 
errors  in  the  pointing  of  the  detector,  sufficient  care  was  taken 
in  obtaining  the  measurements  so  as  to  avoid  these  particular 
problems.  The  measurements  were  found  to  be  very  repeatable  so 
as  to  rule  out  changing  flashlamp  intensities  and  scintillation 
errors.  Although  the  attenuation  coefficient  may  decrease  with 
distance  due  to  preferential  attenuation  of  the  lower  wavelength 
light  within  the  band,  or  due  to  spatial  nonuniformities  of  the 
atmosphere,  the  observed  inconsistencies  were  opposite  to  those 
expected  from  these  effects.  The  major  contributor  to  the  incon- 
sistencies was  found  to  be  the  changing  visibility  during  the  time 
the  measurements  were  obtained. 

Several  steps  were  taken  to  remedy  this  problem.  Measurements 
were  subsequently  made  in  the  afternoon  as  opposed  to  before  noon, 
when  the  early  morning  haze  was  still  being  "burned  off."  The 
detector  pointing  process  was  greatly  simplified  to  save  time. 
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A hood  was  placed  over  the  lamp  to  eliminate  scattered  radiation 
in  all  but  the  immediate  forward  direction  (see  Figure  21),  Attenua- 
tion measurements  were  completed  as  quickly  as  possible  before  any 
scattering  measurements  were  begun.  In  this  manner,  completely 
consistent  measurements  were  made  between  various  sites. 

Experiments  have  been  carried  out  to  simultaneously  measure 
the  attenuation  coefficient  at  three  different  wavelengths  using 
the  filters  described  In  Figure  13.  The  attenuation  coefficients 
measured  using  the  2680,  2580,  and  2^80  filter  combinations  were 
2.65,  3-05*  and  3.05  respectively.  Visibility  was  observed  to 
be  between  5 and  10  km,  which  is  consistent  with  the  data  shown 
In  Figure  22.  It  Is  worth  noting  In  Figure  13  that  the  three 
filters,  when  used  In  combination  with  filters  no.  1 and  2,  do 
not  peak  at  the  same  wavelength  as  they  do  in  Figure  12. 

Scattering  Measurements 

The  procedure  for  mapping  atmospheric  scattering  footprints 
Is  shown  schematically  In  Figure  23.  The  source  was  located  at  some 
slightly  elevated  position,  such  as  the  roof  of  a relatively  low 
building,  and  the  detector  was  moved  some  distance  x away.  In  all 
cases,  the  detector-to-source  line  of  sight  was  within  3*5  degrees 
of  the  horizontal  and  was  defined  as  zero  azimuth  and  zero  elevation. 
Typically,  the  detector  optical  axis  was  raised  vertically  to  some 
angle  of  elevation  denoted  as  ELEV  and  then  rotated  In  the  horizontal 
plane  to  some  azimuthal  position  denoted  as  AZIM.  The  flashlamp  was 
flashed  from  3 to  5 times  and  the  resulting  detector  photon  counts 
were  recorded.  The  procedure  was  repeated  with  the  detector  optical 
axis  systematically  rotated  to  various  azimuthal  angles  and  raised 
to  various  elevations. 

For  each  position  of  the  detector,  optical  axis  background  was 
measured  by  observing  the  photon  count  level  between  flashes  of  the 
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Figure  23.  Geometry  of  atmospheric  scattering 
measurements. 
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lamp.  The  data  were  later  reduced  by  subtracting  the  averaged 
background  count  from  the  average  flashlamp  count  at  each  position. 
On  some  occasions,  the  detector  was  moved  to  a second  location  and 
the  entire  procedure  was  repeated. 

The  very  first  series  of  measurements  showed  that  the  atmos- 
pheric scattering  footprints  were  symmetric  with  respect  to  positive 
and  negative  azimuthal  angles.  In  other  words,  the  averaged  re- 
duced photon  count  at  (20°  ELEV,  AZIM)  was  equal  to  the  count 

at  (20°  ELEV,  -il5°  AZIM)  and  so  forth.  Subsequent  measurements 
were  therefore  limited  to  azimuthal  angles  between  either  0 and 
180°  or  0 and  -180°  but  not  both. 

Samples  of  atmospheric  scattering  data  obtained  on  two  separate 
occasions  at  a distance  x = .82  km  are  shown  In  Figure  2k.  Data 
were  not  obtained  for  zero  elevation  because  of  nearby  bushes, 
trees,  and  other  natural  or  manmade  obstacles.  The  difference  In 
magnitude  of  the  photon  counts  on  the  two  days  was  due  In  part  to 
a change  in  the  setting  of  the  pulse  amplitude  discriminator,  and 
In  part  to  different  weather  conditions.  If  the  data  of  Figure  2k 
are  plotted  as  a function  of  the  absolute  angle  y > between  the 
detector  optical  axis  and  the  detector  to  source  line  of  sight,  we 
obtain  the  results  shown  In  Figure  25. 

The  symmetry  with  respect  to  y shown  In  Figure  25  was  some- 
what surprising  when  first  observed.  Prior  to  these  experiments. 

It  was  expected  that  the  Intensity  of  scattered  light  would  fall 
off  at  lower  elevations  due  to  absorption  effects  by  the  ground. 

It  was  later  found  that,  due  to  strong  preferential  scattering 
in  the  forward  direction,  most  of  the  photons  reaching  a given 
scattering  volume  are  traveling  very  nearly  In  a radial  direction 
from  the  source.  This  simple  observation  was  easily  confirmed  by 
the  high  photon  counts  observed  at  the  lower  angles  in  Figure  25. 

Due  to  the  symmetry  with  respect  to  the  absolute  angle  y , many 


experiments  were  restricted  to  an  azimuthal  angle  of  zero  with 
elevations  ranging  from  0 to  160°. 

Attempts  were  made  on  several  occasions  to  obtain  data  at 
more  than  ore  source  to  detector  distance  x . Unfortunately, 
much  of  this  data  had  to  be  discarded  due  to  the  inconsistencies 
, In  the  measured  attenuation  coefficients  as  reported  In  the  last 

section.  Since  the  primary  source  of  the  inconsistencies  was  Iden- 
tified as  changing  atmospheric  conditions,  no  meaningful  conclusions 
could  be  drawn  from  the  associated  scattering  data.  In  an  attempt 
to  minimize  the  possibility  of  changing  atmospheric  conditions,  sub- 
f sequent  experiments  were  scheduled  In  the  afternoons  to  avoid  the 

; burning  off  periods  of  early  morning  hazes  and  fogs. 

The  first  two  attempts  to  obtain  scattering  data  In  the  after- 
noon at  two  different  x locations  were  ruined  by  sudden  large  In- 
creases In  the  detector  dark  count.  The  source  of  the  high  dark 
counts  was  found  to  be  the  heat  of  the  afternoon  sun.  The  black 
lens  holders  (Figure  lA)  which  held  the  solar  blind  filters  were 
becoming  hot  to  the  touch  and  this  heat  was  being  transferred  to 
the  photomultiplier  tube.  Once  the  lens  holders  were  painted 
white  and  a small  sun  shade  vras  Installed  over  the  detector  box, 
this  type  of  increase  in  background  count  was  no  longer  observed. 

Finally,  a procedure  for  obtaining  scattering  measurements 
and  the  associated  attenuation  coefficient  was  devised  which  gave 
consistent  results.  First,  scattering  measurements  were  made  at 
the  furthest  location  (2.08  km),  and  then  a source-to-detector  line 
of  sight  reading  was  obtained  using  the  stop  and  lens  as  described 
In  the  apparatus  section.  The  detector  equipment  was  then  quickly 
moved  to  the  nearest  location  ( .82  km)  and  another  line  of  sight 
measurement  obtained.  Finally,  scattering  measurements  were 
obtained  at  the  nearest  location.  The  entire  procedure,  from  the 
first  measurement  at  the  furthest  location  to  the  last  measurement 


44 


KH-. OKT*' mti/^f.'^ yyc --wsti^^pyrs* »■*"'* -in  •,->.< 


5^.  the  nearest  location,  required  approximately  ^0  minutes.  The 
results  of  one  such  measurement  are  shown  In  Figure  26.  The 
associated  attenuation  coefficient  was  found  to  be  2.65/km.  At 
the  2.08  km  Docatlon,  data  were  obtained  at  more  than  one  azi- 
muthal angle  In  order  to  check  the  symmetry  already  observed  at 
the  .82  km  location.  The  symmetry  appeared  to  be  maintained, 
although  the  photon  count  for  angles  greater  than  60°  was  less 
than  10  counts  and  therefore  subject  to  scintillation  error  from 
the  background  count,  which  was  of  the  order  of  10-12  counts/sec. 
The  dashed  and  solid  lines  on  Figure  26  represent  numerical  solu- 
tions to  a simplified  atmospheric  scattering  model  which  will  be 
presented  In  the  next  section. 

Aluminum  Reflector  Experiment 

An  experiment  was  carried  out  to  determine  the  effects  of 
placing  a small  aluminum  reflector  under  the  flashlamps.  Intui- 
tively, It  was  believed  that  this  reflector  would  Increase  the 
observed  photon  count  by  Increasing  the  amount  of  UV  radiation 
entering  the  atmosphere.  At  close  range  (x  » .2  km),  the  photon 
counts  were  observed  to  Increase  significantly.  The  data  obtained 
at  X ■ .82  km,  with  and  without  the  aluminum  reflector,  are  shown 
In  Figure  27.  As  can  be  seen,  the  photon  count  remained  approxi- 
mately the  same,  whether  the  reflector  was  used  or  not.  These 
unexpected  results  were  verified  by  repeating  the  experiment 
several  times  during  a very  short  period  of  clme.  It  was  later 
observed,  from  a series  of  experiments  discussed  below,  that  there 
Is  a strong  favoring  of  scattering  In  the  forward  direction. 
Consequently,  very  little  of  the  radiation  leaving  the  flashlamp 
In  the  vertical  direction  (the  direction  of  light  reflected  by 
the  plate)  will  later  reach  a point  In  the  horizontal  plane  of  the 
lamp.  Conversely,  the  scattered  light  reaching  the  deflector  can 
be  significantly  Increased  by  rotating  the  aluminum  reflector  into 
the  vertical  plane  behind  the  flashlamp,  thereby  reflecting 
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additional  UV  light  Into  the  horizontal  plane  In  the  direction 
of  the  detector. 

Blocking  Experiments 

The  purpose  of  this  series  of  experiments  was  to  determine 
the  contributions  of  various  portions  of  the  atmosphere  to  the 
signal  reaching  the  receiver.  The  forward  preference  of  the 
single  scattering  phase  function  played  an  Integral  part  In  these 
experiments.  Unfortunately,  the  single  scattering  phase  functions 
In  the  UV  have  never  been  measured.  One  of  the  purposes  of  these 
measurements  was  to  verify  the  assumed  mathematical  formulation  of 
the  phase  function,  which  will  be  discussed  In  the  next  section. 

The  experiments  were  carried  out  by  using  a device  which 
Intercepted  all  of  the  emitted  radiation  from  the  flashlamp  up  to 
a given  angle  from  the  horizontal,  as  shown  In  Figure  28.  The 
Interception  angle  varied  from  -3°  to  +22°  In  Increments  of  appro- 
ximately 2.5°.  While  zero  blocking  angle  corresponds  to  the 
source-to-detector  line  of  sight,  negative  blocking  angles  were 
used  because  of  the  finite  thickness  of  the  flashlamp.  For  each 
blocking  angle,  the  detector  angle  was  varied  from  0 to  l60°. 

The  overall  results  of  one  such  survey  are  shown  In  Figure  29.  It 
can  be  seen  from  Figure  29  that,  while  the  average  drop  in  photon 
count  Is  nearly  50!?  for  the  maximum  blocking  angle,  the  decrease 
Is  actually  less  than  the  decrease  associated  with  rotating  the 
dectector  an  equal  number  of  degrees  further  away  from  the  source. 
The  50%  overall  decrease  Is  also  equivalent  to  moving  the 
detector  an  additional  .25  km  away. 

Ideally,  It  would  have  been  desirable  to  perform  the  blocking 
experiments  all  the  way  up  to  a blocking  angle  of  90°,  but  a 
considerably  larger  and  more  complicated  blocking  device  would 
have  been  required.  Nevertheless,  a sufficient  trend  can  be  seen 
from  Figure  29  to  conclude  that  transmitting  and  receiving  from 
foxholes  or  from  opposite  sides  of  a hill  Is  certainly  possible. 
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Figure  28.  Geometry  of  blocking  experiments. 


Effects  of  Rain  and  Scintillation 


Due  to  a fortunate  turn  of  events,  it  was  possible  to  deter- 
mine the  effects  of  rain  on  the  atmospheric  scattering.  Data 
were  obtained  for  two  consecutive  days  under  virtually  identical 
conditions  including  location,  lamp  intensity,  voltage  settings, 
visibility,  temperature,  time  of  day,  etc.  During  the  second 
day,  a light  misty  rain  began  to  fall  such  that  the  visibility 
was  reduced  to  somewhere  between  1.5  and  2.0  km.  The  measured 
photon  counts  increased  dramatically,  as  shown  in  Figure  30.  The 
intensity  of  the  rainfall  was  such  that  it  took  approximately 
five  minutes  for  a data  sheet  to  become  too  wet  for  writing 
purposes.  The  increase  in  signal  is  significant  because  it  shows 
that  communications  using  ultraviolet  scattering  would  improve 
during  bad  weather  while  performance  of  other  systems  might  be 
degraded. 

Scintillation  measurements  were  made  on  three  different 
occasions,  with  the  results  of  two  shown  in  Figure  31  and  Table  2. 
It  has  been  found,  however,  that  a good  portion  of  the  apparent 
scintillation  shown  in  Figure  31  is  due  to  electronic  noise 
resulting  from  Inadequate  shielding  of  the  detector  electronics. 

It  can  be  seen  from  Table  2 that  scintillation  effects  Increase 
with  distance  as  might  be  expected. 


51 


^ -tn*  •y^ 


i 

I 

I 

1 

[ 

t 

[ 

[ 


IV.  THEORETICAL  INVESTIGATION 

The  analysis  of  atmospheric  scattering  of  UV  radiation 
can  be  approached  in  a number  of  ways.  Case  and  Zweifel  (Ref. 

4)  have  investigated  the  isotropic  scattering  process  In  an 
infinite  atmosphere.  Levine  (Ref.  5)  analyzed  Isotropic  scatter- 
ing in  a uniform  atmosphere  of  finite  depth.  Initially,  our 
efforts  v/ere  devoted  to  solving  the  radiative  transfer  equation 
for  the  more  general  case  of  isotropic  scattering  with  absorption 
and  with  a finite  field  of  view  detector.  The  results  of  this 
initial  effort  are  given  in  Reference  6.  Because  of  certain 
limiting  assumptions,  such  as  isotropic  scattering  and  weak 
absorption,  the  solution  to  the  radiative  transfer  equation 
given  in  Reference  6 does  not  agree  with  observed  data. 

Simplified  Scattering  Model 

The  solution  of  the  radiative  transfer  equation  with  the 
additional  complication  of  a nonlsotropic  scattering  function  is 
beyond  the  scope  of  the  present  research.  An  attempt  has  been 
made  to  construct  a simplified  seml-emplrlcal  model  which  will 
accurately  predict  the  observed  atmospheric  scattering  footprints. 
The  basic  approach  in  the  development  of  this  seml-emplrlcal  model 
was  to  account  for  as  many  intuitive  factors  and  observed  physical 
phenomena  as  possible.  Only  a few  constants  have  been  introduced 
into  the  model,  all  of  which  have  some  physical  significance. 

The  criteria  which  have  been  established  for  developing  the 
simplified  model  are  given  below.  Basically,  the  model  should 

1.  Account  for  atmospheric  absorption  of  the  UV 
radiation  (primarily  due  to  ozone). 

2.  Account  for  absorption  by  the  earth's  surface. 

3.  Give  higher  photon  counts  at  close  range  as  the 
scattering  coefficient  is  increased.  This  phenomenon 
was  observed  as  shown  in  Figure  30. 
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Give  lower  photon  counts  at  long  range  as  the 
scattering  coefficient  is  Increased.  Although 
this  phenomenon  has  not  actually  been  observed, 
it  makes  sense  intuitively  that  with  increased 
scattering  the  total  path  length  traveled  by  a 
photon  in  reaching  a location  far  from  the  source 
must  increase.  Since  ozone  absorption  effects  are 
proportional  to  total  path  length,  the  net  effect 
would  seem  to  be  a decrease  in  the  number  of  photons 
reaching  the  far  location. 

5.  Utilize  a generalized  scattering  function  with 
arbitrary  constants  which  control  the  preference 
for  scattering  In  the  forward  direction.  The 
surface  integral  of  the  function  over  steradlans 
should  be  equal  to  one. 

6.  Give  the  correct  shape  of  the  atmospheric  scattering 
footprints  observed  in  the  field. 

7.  Give  the  correct  magnitude  of  the  scattering  foot- 
prints simultaneously  at  various  x locations  (see 
Fig.  23). 

8.  Give  the  correct  magnitude  of  photon  counts  when  the 
source  lamp  Is  within  the  field  of  view, 

9.  Give  the  correct  Intensities  for  the  blocking  experiments. 

10.  Incorporate  the  measured  fleld-of-view  function. 

Assumptions  and  Model  Development 

Most  of  the  UV  photons  reaching  a point  in  the  atmosphere  are 
traveling  very  nearly  In  a radial  direction  from  the  flashlamp. 

This  assumption  Is  easily  verified  by  looking  at  Figure  26  which 
shows  the  angular  dependence  of  the  radiation  reaching  a point  In 
space  occupied  by  the  detector  Itself.  Two  phenomena  are  believed 
to  be  responsible  for  this  effect.  The  basic  single  scattering 
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phase  function  is  strongly  favored  in  the  forward  direction. 
Also,  photons  which  travel  more  indirect  routes  are  more  likely 
to  be  absorbed  by  atmospheric  ozone. 


The  basic  geometry  used  in  the  model  is  given  in  Figure  23. 
Each  point  in  space  is  defined  by  a vector  ^ which  originates 
from  the  flashlamp.  Each  point  in  space  is  also  defined  by  the 
vector  r which  originates  from  the  detector.  In  the  absence 
of  scattering,  the  photons  which  would  pass  through  a small  area 
AA  , which  is  perpendicular  to  R , would  be  equal  to 


I AAe 
o 


-°ABS” 
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where  is  the  total  number  of  photons  per  flash  radiated 

isotropically  into  space,  is  the  absorption  coefficient, 

and  R is  the  length  of  vector  R . The  “cerm  AA/MttR^  Is,  of 
course,  the  steradlancy  of  the  area  AA  with  respect  to  the 
flashlamp . 

The  fraction  of  the  photons  which  are  scattered  before 
reacing  the  location  at  R is  given  by 

“°SCATT^ 

e 

where  ^SCATT  scattering  coefficient.  Of  those  which  are 

scattered,  n*  are  scattered  into  the  backward  direction  and 
1 - n'  are  scattered  into  the  forward  direction.  It  is  assumed 
that  those  scattered  in  the  forward  direction  will  still  reach 
the  location  R . Thus,  the  total  number  of  scattered  and  trans- 
mltted  photons  reaching  the  location  R is  given  by 

^ 
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The  fraction  of  photons  n*  scattered  Into  the  backward  direction 
will  be  determined  later  by  integrating  the  single  scattering 
phase  function  over  the  rearward  facing  hemisphere. 


Of  those  photons  passing  through  the  area  AA  , the  fraction 
scattered  over  a length  AH  (which  is  parallel  to  R)  Is,  by 
definition,  ‘^SCATT  ’ those  which  are  scattered,  the 

fraction  scattered  per  unit  solid  angle  in  the  direction  of  the 
detector  is  given  by  the  single  scattering  phase  function  P(0,a) 
where  3 is  the  angle  the  vector  r makes  with  the  vector  R as 
shown  in  Figure  23,  and  ”a"  is  a parameter  which  represents  the 
scattering  characteristics  of  the  existing  aerosol  distribution. 
Thus,  the  total  number  of  photons  scattered  in  the  volume  element 
AA  • AS,  and  heading  In  the  direction  of  the  detector  is  given  by 


^Trr‘ 
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where  Aj^/^irr  is  the  steradlancy  of  a detector  of  area  Aj^ 
with  respect  to  the  scattering  volume  which  is  located  at  a dis- 
tance r (see  Fig.  23). 


Not  all  of  the  photons  leaving  the  scattering  volume  will 
arrive  at  the  detector,  due  to  absorption  and  scattering.  Using 
the  same  reasoning  as  was  applied  for  scattering  and  absorption 
along  R , the  fraction  of  photons  lost  along  r would  be 

■■^(‘^abs'‘'’^’®scatt) 
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Of  the  photons  reaching  the  detector,  only  a fraction  will  actually 
result  in  photon  counts  due  to  the  effects  of  filter  transmission, 
photomultiplier  tube  efficiency  and  geometric  losses.  All  of  these 
effects  are  accounted  for  by  the  fleld-of-vlew  function  P0V(6) 
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where  0 is  the  angle  between  the  detector  optical  axis  and 
the  vector  r In  Figure  23.  The  angular  dependence  of  the  fleld- 
of-vlew  function  for  the  detector  used  In  this  study  has  been 
measured  and  is  shown  In  Figure  15.  The  total  photon  count  PC 
due  to  all  volume  elements  AV  = AA  • Mi  In  the  fleld-of-vlew 
cone  of  the  detector  Is  then  given  by 


PC(y) 


■/// 


"(^■’■r)  (°^ABS‘*t^  *°SCATt) 


TT'srFT 

Idtt  R r 


AV  (5) 


Equation  (5)  is  the  basic  relation  used  for  the  simplified 
scattering  model.  Because  of  its  complicated  form,  Eq.  (5)  has 
been  Integrated  numerically  by  computer  program.  Some  of  the 
algebraic  relations  necessary  to  Integrate  the  equation  are  given 
later  In  this  section. 


Single  Scattering  Phase  Function 

According  to  the  results  of  a fairly  extensive  literature 
survey,  single  scattering  phase  functions  In  the  2200  to  2800X 
region  of  the  spectrum  have  never  been  measured.  We  therefore 
initially  chose  the  Henyey-Greensteln  phase  function  given  by 


P(B,g) 


1 


(1  + 


- 2 cos 


(6) 


where  g is  the  average  value  of  the  cosine  of  the  scattering 
angle.  The  Henyey-Greensteln  function  represents  a general  family 
of  curves  with  a forward  scattering  preference  as  shown  In  Figure 
32.  (Note  that  the  n In  Figure  32  Is  equal  to  1 - n’  .)  These 
functions  were  chosen  because  they  exhibit  some  of  the  qualitative 
features  of  the  more  realistic  Delrmendjian  phase  functions  shown 
in  Figure  33.  The  Deirmendjlan  functions  were  obtained  by 
averaging  the  Mle  scattering  functions  over  polydispersed  aero- 
sols. The  polydispersed  aerosols  are  classified  as  haze  types. 
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Figure  33.  DelrmendJ Ian '3  phase  functions. 
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as  shown  in  Piguro  3^  where  number  density  is  plotted  as  a 
function  of  particle  radius.  The  Delrmendjlan  functions  shown 
in  Figure  33  correspond  to  haze  type  L of  Figure  3^.  The 
primary  qualitative  difference  between  the  Henyey-Greensteln  and 
Deirmendjian  phase  functions  is  in  the  region  of  scattering 
angles  beyond  60°  as  shown  in  Figure  35. 

Because  the  Delrmendjlan  functions  exist  only  in  the  form  of 
tabulated  values  for  specific  wavelengths  and  haze  types,  and 
because  they  have  not  been  tabulated  for  wavelengths  below  ^500&, 
an  attempt  has  been  made  to  modify  the  Henyey-Greenstein  functions 
to  look  more  like  Delrmendjlan  functions.  The  modified  Henyey- 
Greensteln  functions,  which  for  simplicity  will  be  referred  to  as 
the  Neer-Sandrl  functions,  are  written  as 


p(e,g) 


1 

w 


g ~-S  ^ + X'(3  C0S^3-1) 

L(1  + g*^  - 2 cos 


(7) 


where  the  second  term  is  the  second  Legendre  polynomial.  Because 
the  second  term  is  symmetric  about  B " 90° , the  Integral  of 
P(B,g)  over  all  solid  angles  is  unchanged.  In  the  same  fashion, 
g , the  average  value  of  cos  6 , is  also  unchanged. 


Since  the  phase  function  must  alvays  be  greater  than  or 
equal  to  zero,  a value  of  X'  must  be  chosen  which  will  ensure 
that  the  minimum  value  of  P(B,g)  is  positive.  By  setting  tiie 
derivative  of  Eq.  (7)  equal  to  zero,  we  can  express  X in  terms 
of  the  scattering  angle  B^  which  corresponds  to  the  minimum 
value  of  P(B»[  ). 


X 




2 cos  BgCl  + g^-S  cos  Bq) 


(8) 


Equation  (7)  can  then  be  written  as 
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xa  inr;  depend  on  the  particular  mod#»i 
(Reproduced  from  Reference  8) 
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p(^) 


P(3,g) 
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1 - 


(1  + g - 2gcos  e) 
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g(l  - g^)(3  cos^e  - 1) 

2 1 cos  0^1  ( 1 + g^  + 2g  I c os  3q  I ) 


(9) 


A value  of  1/7  has  been  chosen  for  (cos  6^)  because  it  not  only 
guarantees  that  the  phase  function  will  be  everywhere  positive 
but  also  causes  the  N'eer-Sandri  function  to  look  like  Diermend- 
Jlan's  function  as  shown  in  Figure  35* 


Multiple  Scattering  Effects 

The  atmospheric  scattering  model  described  by  Eqs . (5)  and 
(9)  has  one  obvious  fault.  As  the  distance  R increases,  the 
number  of  times  that  the  avei’age  photon  has  been  scattered  before 
reaching  that  distance  Increases.  The  assumption  that  all  photons 
are  traveling  parallel  to  R becomes  less  valid  as  R Increases. 
Thus,  the  direction  of  those  photons  scattered  at  R will  not 
favor  the  forward  direction  as  much.  In  other  words,  the  net 
effect  is  that  the  scattering  would  appear  to  become  more  iso- 
tropic as  the  distance  R Increases.  One  way  to  account  for  this 
effect  is  to  allow  the  effective  g value  (where  g is  the 
average  value  of  the  cosine  of  the  scattering  angle)  to  decrease. 
The  effect  of  decreasing  g can  be  seen  in  Figure  32  where,  as 
g decreases,  the  scattering  function  becomes  more  isotropic.  In 
order  to  account  for  the  effects  of  multiple  scattering,  the 
following  relation  for  an  effective  g value  has  been  assumed: 


g 


_ ^"^l^SCATT^ 


(10) 


where  g^  Is  the  actual  value  for  the  single  scattering  phase 
function  and  c^  Is  a constant  describing  how  fast  the  effective 
g value  should  drop.  The  functional  form  of  Eq.  (10)  has  been 
assumed,  because  for  low  values  of  R most  of  the  photons  will 


not  have  undergone  any  scattering  prior  to  reaching  R , and 
hence  g should  equal  g^  . For  very  large  R , g will  go 
to  zero,  thus  representing  totally  Isotropic  radiation.  It 
would  not  be  unreasonable  to  expect  totally  diffuse  radiation 
for  large  values  of  R , If  It  were  not  for  the  exponential  decay 
with  total  path  length  due  to  ozone  absorption. 


Since  Eq.  (10)  has  been  written  almost  arbitrarily,  it  is 
futile  to  argue  the  functional  dependencies  of  c^  . If,  however 
c^  could  be  expressed  as  a function  of  g^  , and  g^  could  In 
turn  be  measured,  Eq.  (10)  would  Introduce  no  arbitrary  constants 
As  It  stands  now,  c^^  and  g^  are  both  arbitral y constants. 

Now  that  the  scattering  function  has  been  defined  by  Eq.  (9) 
It  Is  possible  to  calculate  n’  which  appears  In  Eq.  (5).  n’ 

Is  defined  as  that  fraction  of  the  scattered  light  which  Is 
scattered  Into  the  rearward  direction.  Hence, 
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According  to  Eq.  (10),  the  effective  value  of  g varies  with  R 
Therefore,  at  each  position  R we  define  a value  of  n*  to 
represent  the  average  value  of  backscattering  along  the  path  from 
0 to  R , according  to  the  relation 


n’(R)  - 


n ' (g«)  + n ' (g(R) ) 


(12) 


Integration  of  Atmospheric  Scattering  Relation 

The  atmospheric  scattering  model  defined  by  Eqs.  (5)  and  (9) 
Is  nearly  Imposslole  to  Integrate  In  closed  form.  Consequently, 
the  equation  has  been  Integrated  numerically  by  dividing  the 


fleld-of-vlew  cone  Into  a finite  number  of  volume  elements  and 
summing.  When  quick  calculations  are  desired,  the  viewing  cone 
can  be  divided  into  circular  slices  which  are  perpendicular  to 
the  optical  axis,  as  shown  in  Figure  36.  The  volume  of  each  slice 
can  be  written  as 

“ Tr(r^  tan  0)^Arjj^  (13) 

and  the  distance  K can  be  expressed  in  terms  of  r,  AZIM,  and 
ELEV  (see  Pig.  23)  as 

R **  J - 2xr  cos  y + (1^) 

The  value  of  cos  6 in  Figure  23  is  found  to  be 

003  6 - rl  - r (15, 


where 

cos  Y - cos(ELEV)  cos(AZIM)  (I6) 

I 

If  more  accuracy  is  desired,  the  fleld-of-vlew  cone  is 
divided  into  elpherlcal  volume  elements,  as  shown  in  Figure  37. 

The  volume  of  these  elements  is  given  by 

- r^  sin  6 Ar  A0  A<|)  (17) 

where  the  more  general  coordinate  system  shown  in  Figure  38  must 
be  used  to  define  6 and  . As  can  be  seen  from  the  figure, 
the  2 coordinate  is  defined  as  the  detector  line  of  sight.  The 
flashlarap  or  source  is  located  a distance  from  the  detector 

at  coordinates  absolute  angle  y away  from 

the  detector  line  of  sight.  Just  as  in  Figure  23.  The  centers 
of  the  spherical  volume  elements  are  defined  in  spherical  coordi- 
nates by  r,  G , and  . The  distance  R from  the  detector  to 
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Figure  37-  Spherical  volume  elements. 
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Figure  38.  General  coordinate  system. 


a specific  spherical  volume  element  is  given  by 


/^2  — 2 

- 2x^r[sln  y sin  0 cos(a-<j!)  + cos  y cos  6]  + r (18) 

and 

X [sin  y sin  0 cos(a-4»)  + cos  y cos  0]  - r 
cos  e •=  ^ (19) 

The  height  above  the  ground  of  a volume  element  can  be  found  If 
It  is  assumed  that  the  detector,  the  lamp,  and  the  y-axls  of  the 
detector  all  lie  within  the  ground  plane.  In  this  event,  the 
height  z is  given  by 

o 

z - r cos  0 sln(ELEV)  - r sin  9 cos  (j)  cos(ELEV)  (20) 

O 

The  computer  program,  which  numerically  Integrates  Eq.  (5), 
calculates  z^  for  each  volume  element  and  registers  a zero  If 
z Is  negative.  When  the  lamp  Itself  1s  within  the  field  of  view, 

o 

the  size  of  the  volume  elements  greatly  decreases  In  the  vicinity 
of  the  lamp.  Singular  point  problems  are  minimized  by  making  sure 

O 

that  the  quantity  (4irR  )/3  is  always  bigger  than  the  volume 
element  Itself. 

Regardless  of  whether  the  cone  slices  of  the  spherical  volume 
elements  are  used,  the  computer  program  Input  includes  '^scatt  * 
^ABS  * , X , ELEV  , AXIM  , and  . For  comparison  with 

atmospheric  scattering  data,  the  value  of  the  attenuation  coeffi- 
cient Is  measured.  In  principle,  should  also  be  determined 

Independently  by  measuring  the  ozone  concentration.  In  this  study, 
a value  of  “ ‘T/km  has  been  assumed,  based  on  tabulated 

values  from  Reference  7.  The  scattering  coefficient  la  then 
determined  from  Eq.  (4).  The  number  of  arbitrary  constants  Is 
thus  limited  to  g^  and  c^  . The  results  of  the  simplified 
analyses  for  “ .5  and  c^  ■ .3  are  compared  with  scattering 
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data  in  Figure  26.  Results  are  also  shown  in  Figure  26  for  the 
case  where  the  Henyey-Greenstein  function  (Eq.  (6))  is  used 
instead  of  the  Neer-Sandri  phase  function  (Eq.  (9)).  The  g^  “ 

.5  and  c^  « .3  case  is  also  compared  to  the  "blocking" 
experiments  in  Figure  29.  The  computer  program  can  be  modified 
to  approximate  the  effects  of  the  blocking  experiment  for 
detector  angles  (see  Fig.  28)  between  40°  and  l40°. 

The  apparent  agreement  of  the  simplified  model  with  the  field 
measurements  is  encouraging.  As  mentioned  earlier,  it  is  possible 
that  some  cancellation  of  error  may  be  Involved.  However,  the 
agreement  between  the  results  predicted  by  the  model  and  the  data 
of  the  blocking  experiments  indicates  that  we  have  at  least 
succeeded  in  predicting  which  portion  of  the  signal  is  coming  from 
a particular  portion  of  the  f leld-of-vlew  cone.  Also,  the  agree- 
ment obtained  at  2.08  km  for  an  angle  of  0°  (see  Pig.  26)  1s 
equally  encouraging.  At  0°,  the  flashlamp  Itself  is  In  the  field 
of  view,  and  the  photon  count  la  calculated  from  the  sum  of  Eqs. 
(2)  and  (5).  If  the  relative  magnitude  of  scattered  and  trans- 
mitted light  calculated  was  in  error,  the  agreement  observed  in 
Figure  26  could  not  be  obtained. 
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NOTE:  CHAPTERS  V AND  VI  HAVE  BEEN  CLASSIFIED 

CONFIDENTIAL  AND  APPEAR  SEPARATELY  IN 
VOLUME  II  OF  THIS  REPORT. 
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VII.  SUMMARY 


Voice  communication  has  been  accomplished  via  scattered 
ultraviolet  radiation  In  the  solar  blind  region  of  the  spectrum. 

The  transmitter  Is  a frequency  modulated,  high  repetition  rate, 
hydrogen-xenon  UV  flashlamp  operating  at  a carrier  frequency  of 
il8l0  pulses  per  second.  The  receiver  Is  composed  of  solar  blind 
filters  and  a rubidium  tellurlde  photomultiplier  tube.  Good 
quality  voice  communications  have  been  transmitted  around  the 
corner  of  a building  due  to  pure  UV  backscatterlng  from  atmos- 
pheric aerosols.  An  analysis  has  been  performed  which  Indicates 
that,  with  optimization  of  the  various  components,  a back  pack 
field  unit  could  send  and  receive  voice  communication  over  reasonable 
distances,  which  depend  on  the  wavelength,  when  operating  In  a 
purely  multidirectional,  non-llne-of-slght  mode.  Naturally,  heavier 
systems  would  have  a much  longer  range.  The  Increase  In  range 
which  can  be  accomplished  with  directional  or  beamed  transmitted 
UV  Is  not  known.  A statistical  analysis  has  shown  that  with  a 
slight  modification  of  the  electronic  circuitry  random  background 
signals  due  to  fires,  flares,  and  other  ultraviolet  radiation 
sources  can  be  rejected.  The  communications  system  appears  to 
operate  better  In  bad  weather  due  to  higher  concentrations  of 
atmospheric  aerosols. 

A high-energy  flashlamp  and  the  ultraviolet  detector  have 
been  used  to  systematically  measure  the  overall  effects  of 
scattering  and  absorption  In  the  lower  atmosphere  for  light 

o o 

between  2500A  and  2700A  In  the  solar  blind  region  of  the  spectrum. 

The  line  of  sight  attenuation  coef i Iclents , measured  over  the 
.82  and  2.08  km  distances  between  the  ground-based  flashlamp  and 
detector  locations,  were  found  to  be  consistent  with  previously 
published  values.  The  atmospheric  scattering  footprints  obtained 
by  surveying  the  sky  were  observed  to  be  symmetric  with  respect  to 
the  angle  between  the  detector  optical  axis  and  the  detector-to- 
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flashlamp  line  of  sl^ht.  It  has  been  concluded  that  due  to 
both  a strong  preference  for  forward  scattering  and  to  an  expon- 
ential decay  with  pathlength  associated  with  ozone  absorption 
most  of  the  photons  reaching  a given  point  in  the  atmosphere  are 
traveling  very  nearly  in  a radial  direction  from  the  flashlamp. 

A light  misty  rain  was  observed  to  Increase  the  detected  radia- 
tion by  a factor  of  four. 

A simplified  theory  based  on  both  experimental  observation 
and  physical  intuition  has  been  developed  to  model  ultraviolet 
scattering  and  absorption  in  the  lower  atmosphere.  The  model 
utilizes  the  generalized  Henyey-Greensteln  phase  function, 
modified  to  more  closely  resemble  DelrmendJ ian’ s phase  functions 
for  polydlspersed  aerosols.  Multiple  scattering  effects  are 
modeled  by  assuming  the  phase  function  to  become  Isotropic  at 
large  distances  from  the  light  source.  Two  arbitrary  constants, 
introduced  to  represent  the  unknown  ultraviolet  phase  function, 
are  chosen  to  obtain  agreement  with  the  magnitude  and  angular 
dependence  of  measured  radiation  at  a given  source-to-detector 
distance.  Simultaneous  agreement  between  theory  and  experiment 
for  another  source-to-detector  distance,  as  well  as  for  other 
related  experiments,  tends  to  validate  the  choice  of  constants 
as  well  as  the  overall  model. 
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